Theory of water and charged liquid bridges 
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The phenomena of liquid bridge formation due to an applied electric field is investigated. A new 
solution of a charged catenary is presented which allows to determine the static and dynamical 
stability conditions where charged liquid bridges are possible. The creeping height, the bridge 
radius and length as well as the shape of the bridge is calculated showing an asymmetric profile in 
agreement with observations. The flow profile is calculated from the Navier Stokes equation leading 
to a mean velocity which combines charge transport with neutral mass flow and which describes 
recent experiments on water bridges. 
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I. INTRODUCTION 



A. The phenomenon 



The formation of a water bridge between two beakers 
under high voltage is a phenomenon known since over 100 
years [lj . When two vessels brought in close contact and 
a high electric field is applied between the vessels, the 
water starts creeping up the beakers and forms a bridge 
which is maintained over a certain distance as schemat- 
ically illustrated in figure [TJ Due to the voltage applied 
by the vessels the electric field is longitudinally oriented 
inside the cylindrical bridge. It has remained attractive 
to current experimental activities @, Q ■ On the one side 
the properties of water are such complex that a complete 
microscopic theory of this effect is still lacking. On the 
other side the formation of water bridges on nanoscales 
are of interest both for fundamental understanding of 
electrohydrodynamics and for applications ranging from 
atomic force microscopy 4] to electrowetting problems 
Q. Microscopically the nanoscale wetting is important 
to confine chemical reactions &] which reveals an interest- 
ing interplay between field-induced polarization, surface 
tension, and condensation [3, Q . 

Molecular dynamical simulations have been performed 
in order to explore the mechanism of water bridges at 




FIG. 1: The schematic picture of water bridge between two 
beakers. 



the molecular level leading to the formation of aligned 
dipolar filaments between the boundaries of nanoscale 
confinements [9|. A competition was found of orienta- 
tion of molecular dipoles and the electric field leading 
to a threshold where the rise of a pillar overcomes the 
surface tension [8J. In this respect the understanding of 
the microscopic structure isessential to explain such phe- 
nomena in micro-fluidics [10]. The problem is connected 
with the dynamics of charged liquids which is important 
for capillary jets [llf. current applications in ink printers 
and electrosprays [1J, [ljj. Consequently the nonlinear 
dynamics of breakup of free surfaces and flows has been 
studied intensively [lj, [ljf . 

Much physical insight can be gained already on the 
macroscopic scale, where the phenomena of liquid bridg- 
ing is not restricted to water but can be observed in 
other liquids too [16| which shows that it has its origin 
in electrohydrodynamics [171 ] rather then in molecular- 
specific structures. The traditional treatment is based 
on the Maxwell pressure tensor where the electric field 
effects comes from the ponderomotoric forces and due 
to boundary conditions of electrodynamics [18j . This 
is based exclusively on the fact that bulk-charge states 
decay on a time scale of the dielectric constant divided 
by the conductivity, eeo/c, which takes for pure water 
0.14ms. This decay-time of bulk charges follows from 
the continuity of charge density p c = — V • j combined 
with Ohm's law j = erE = — crV</> where the source of the 
electric field is given by the potential V 2 </> = — p c /eeo- 
An overview about the different forces occurring in mi- 
croelectrode structures are discussed in |19| . 

This simple Ohm picture leads to a problem in par- 
tially charged liquids. Following the Ohm picture one 
has a constant velocity or current density of charged par- 
ticles caused by the external field and limited by friction. 
Contrary, for incompressible fluids the total mass flux 
cannot be constant but is dependent on the area where 
it is forced to flow through. Both pictures seem to be im- 
possible to reconcile. Here in this paper we will present 
a discussion of this seemingly contradiction leading to a 
dynamical stability criterion for the water bridge and a 



combined flow expression. This is in line with the idea 
of [l7[ where the bulk charges have been assumed to be 
realized in a surface sheet. While there the migration 
of charges to the surface has been considered forming a 
charged surface sheet, we adopt here the view point of ho- 
mogeneously distributed bulk charges which flow in field 
direction rather than forming a surface sheet. 

In the absence of bulk charges the forces on the water 
stream are caused by the pressure due to the polarizabil- 
ity of water described by the high dielectric susceptibil- 
ity e. This pressure leads to the catenary form of water 
bridge like a hanging chain [20] - While already the sim- 
plified model of [16| employing a capacitor picture leads 
to a critical field strength for the formation of the water 
bridge, the catenary model [2(| has not been reported 
to yield such a critical field. In this paper we will show 
that even the uncharged catenary provides indeed a min- 
imal critical field strength for the water bridge formation 
in dependence on the length of the bridge. This criti- 
cal field strength is modified if charges are present in the 
bridge which we will discuss here with the help of a new 
charged catenary solution. This allows us to explain the 
asymmetry found in the bridge profile [3j. 



B. Overview about the paper 

The scenario of water or other dielectric bridges is 
thought as follows. Applying an electric field parallel 
to two attached vessels the water creeps up the beaker 
and form a bridge as it is nicely observed and pictured 
in 0. This bridge can be elongated up to a critical field 
strength and it forms a catenary which becomes asym- 
metric for higher gravitation to electric field ratios Q- 
The critical value for stability is sensitively dependent 
on ion concentrations breaking off already at very low 
concentrations. The amount of mass flow through the 
bridge does not follow simple Ohmic transport as we will 
see in this paper. The schematic picture of the water 
bridge is given in figure [T] 

In this paper we want to advocate the following pic- 
ture. Imaging a snapshot of the charges flowing through 
the bridge we cannot decide whether the observed charges 
are due to static bulk charges or due to the floating mo- 
tion of Ohmic bulk charges. This flow of charges within 
the liquid bridge we can associate with a dynamical bulk 
charge in the mass motion which is not covered by the 
decay of Ohmic bulk charges discussed above. Such a 
picture is supported by the experimental observation of 
possible copper ion motion [2l| and by the observation 
that the water bridge is highly sensitive to additional ex- 
ternal electric fields 22] . Strong fields even create small 
cone jets [2J. This dynamical bulk charge will lead us 
to the necessity to solve the catenary problem including 
bulk charges. Though c harg ed membranes have been dis- 
cussed in the literature [23j , a new analytical solution of 
the charged catenary is discussed in this paper. 

The picture of Ohmic resistors and capacitor as de- 



density p = 10 3 kg/m 3 

dielectric susceptibility e = 81 

surface tension a s = 7.27xlO _2 N/m 

viscosity r\ = 1.5xlO _3 Ns/m 2 

conductivity of 

clean water <ro = 5xlO _6 A/Vm 

molecular conductivity 

of NaCl A = 12.6xKr 3 Am 2 /Vmol 

heat capacity c v = 4.187 J/gK 

TABLE I: Variables and parameters used within this paper 
for water. 



scribed above is not sufficient, as one can see from the 
observation that adding a small amount of electrolytes to 
the clean water destroys the water bridge almost imme- 
diately. In other words good conducting liquids should 
not form a water bridge. We will derive an upper bound 
for charges possibly carried in water in order to remain in 
stable liquid bridges. Though we present all calculations 
for water parameters summarized in table HI the theory 
applies as well to any dielectric liquid in electric fields. 

Four theoretical questions have to be answered: (i) 
How is the electric field influencing the height z max wa- 
ter can creep up? (ii) What is the radius R(x) along the 
bridge? (iii) What is the form z — f(x) of the water 
bridge? What are the static constraints on the bridge? 
(iv) Which dynamical constraints can be found for pos- 
sible bridge formation? 

We will address all four questions with the help of four 
parameters composed of the properties summarized in 
table Q] of water. The first one is the capillary height 



3.8 mm 



(1) 



with the surface tension a s , the particle density p and 
the gravitational acceleration g. The second parameter is 
the water column height balancing the dielectric pressure 
called creeping height in the following 



b(E) = 



pg 



= 7,22E' cm 



(2) 



where the dimensionless electric field E is in units of 
10 4 V/cm. The third one is the dimensionless ratio of 
the force density on the charges by the field to the grav- 
itational force density 



c{Pc,E) 



p c E 

pg 



l5.97Ep c 



(3) 



where the charge density p c is in units of ng/l. For dy- 
namical consideration the characteristic velocity 



u 



pgg 

3277 



3.02m/s 



(4) 



will be useful as the fourth parameter. 

The outline of the paper is as follows. In the next sec- 
tion we repeat shortly the standard treatment of creeping 



height and bubble radius of a liquid but add the pressure 
by the external electric field on the dielectric liquid. Then 
we present the form of the bridge in terms of a new so- 
lution of the catenary equation due to bulk charges in 
section IV. In section V we present the flow calculation 
proposing the picture of moving charged particles due to 
the field which drag the neutral particles. This will lead 
to a dynamical stability criterion. Then we compare with 
the experimental data and show the superiority of the 
present treatment. Summary and conclusion ends up the 
discussion in section VI. 



II. ANSWER TO QUESTION (I): CREEPING 
HEIGHT AND (II): RADIUS OF BRIDGE 

We start to calculate the possible creeping height and 
use the pressure tensor for dielectric media [18| 
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FIG. 2: The schematic picture of water bridge creeping up 
the vessel due to the applied electric field. 
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where p is the pressure in the system, R\ , R2 the principal 
radii of curvature such that the second term on the right 
hand side describe the contribution due to surface tension 
and the last terms are the parts due to the forces in the 
dielectric medium. We assume a density-homogeneous 
liquid such that for the dielectric susceptibility e = e — 
p{de/ dp)T ~ £• Further we consider first the stationary 
problem which means that viscous forces can be neglected 
in ©. 

Denoting the components of the normal vector by e k , 
the stability condition between water (W) and air (A) is 
given by 



^.ik k 
a (A) e (A) 



ik k 
T (W) e (W) 



rr ik p k 

a (A) e (W)- 



(0) 



Since the principal curvature of the tube is much larger 
radially than parallel, we have R2 ~ 00 and denoting 
the coordinate in the direction of the height with z, the 
pressure difference between water and air is p\y — pl = 
pgz. We employ the boundary conditions for the normal 
E n and tangential E f components of the electric field 



E {A) - eE {W) — eE n, ^(A) 

and the balance © with ([5]) reads 



rpt rpt rp 
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pgz 



Ri 
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(7) 



(8) 



Please note that due to the migration of charges to the 
surface one should consider a surface charge here in prin- 
ciple. We adopt thorough the paper the simplified pic- 
ture that the charges remain bulk-like due to the pre- 
ferred motion along the field and no surface charges are 
formed. The influence of such surface charges is con- 
sidered as marginal since the curvature of the bridge is 
minimal leading to preferential tangential components of 
electric fields. 



We assume the electric field in x-direction such that 
E t — —Ecosa, E n = Esuia where z'{x) = tana is the 
increase of the surface line of the water as illustrated in 
figure [U Using the parameters ([T]) and (J2J we obtain 
from the stability condition $E§ the differential equation 



2z 



(1 + Z ' 2 ) 3 / 2 



fo(e-l) 

pg 



{eE 2 n + E 2 )^b (9) 



where we used the approximation of small normal electric 
fields justified if there are no surface charges. This shows 
the modification of the standard treatment of capillary 
height by the applied field condensed on the right hand 
side. The first integral of ([9]) is 



VT 



z' 2 a 2 



(10) 



and we have used the condition that for x — > 00 the sur- 
face is z = z' = 0. The explicit solution of the surface 
curve z{x) is quite lengthy and not necessary here. In- 
stead we can give directly the maximally reachable height 
in dependence on the electric field. Therefore we use the 
angle 6 = 90 — a of the liquid surface with the wall such 
that z'(x) — — cot 6 and from ([TU)) we obtain 



b 2 b 

— + a 2 (l-sm8) < - 
4 2 



a — ^max 



(11) 



which shows that without electric field the maximal 
creeping height is just the capillary length ^ as it is 
well known. The other extreme of very high fields leads 
to the field-dependent length © which justifies the name 
creeping height. This answers the first question concern- 
ing creep heights. 

The second question, how large the radius of the bridge 
is, one finds by equating the pressure due to surface ten- 
sion with the gravitational force density 



R 



pgz sa pg 2R 



(12) 



such that the radius of the water bridge is at the beaker 



B. Static stability criteria 



R « a/2. 



(13) 



Without using this approximation we could express the 
curvature again by differential expressions in z(x) defin- 
ing a radial profile, as it can be found in text books {18j | . 
The radius of the bridge at the beaker is nearly indepen- 
dent on the applied electric field but only dependent on 
the surface tension and gravitational force. Along the 
bridge the radius will change with the applied electric 
field as we will see later in section IV. C. 



III. ANSWER TO QUESTION (III): LIQUID 
BRIDGE SHAPE 



A. Charged catenary 

Now we turn to the question which form the water 
bridge will take. Therefore we consider the center of 
mass line of the bridge being described by z = f(x) 
with the ends at /(0) = f(L) = 0. The force densi- 
ties are multiplied with the area and the length element 
ds = y/l + f' 2 dx to form the free energy. We have the 
gravitational force densitypgf and the volume tension 
pgb as well as the force density by dynamical charges 
p c Ex which contributes. The surface tension is negligi- 
ble here. The form of the bridge will be then determined 
by the extreme value of the free energy 

L L 

?{z)dx = pgfmx) + b - caOVT + T*** - extr. 



(14) 

where c is given by ([3]) and b defined in ((5J). 

As shown in [24j and shortly outlines in appendix [A] 
the solution can be represented parametrically as 



Without dynamical bulk charges, c = 0, d = 1, the 
solution (fT5j) is just the well known catenary [2Q|. The 
boundary condition dl~7l) reads in this case 

lb 2£ L 

— = j- cosh— >£ c = 1.5088... (18) 

which means that without bulk charges the condition for 
a stable bridge is 



b > \l£, c 



(19) 



Together with ([2]) this condition provides a lower bound 
for the electric field in order to enable a bridge of length 
L. This lower bound for an applied field appears obvi- 
ously already for the standard catenary and has been not 
discussed so far. 

Lets now return to the more involved case of bulk 
charges and the new solution of charged catenary (fl~5j) . 
The field-dependent lower bound condition ([T7l) is plot- 
ted in figure [3J One see that in order to complete p7|) 
the bulk charge parameter c has to be lower than the 
maximal value of c m which reads 



c < c m (£ ,b) 



(20) 



and which is plotted in the inset of figure [3j Remember- 
ing the definition of the bulk charge parameter ([3j we 
see that (|2"0|) sets an upper bound for the bulk charge in 
dependence on the electric field. The lower bound ([T9"j) of 
the electric field for the case of no bulk charges is obeyed 
as well since the curve in the inset of figure [3J starts at 
b > L£ c /2 which is the lower bound already present for 
uncharged catenaries (fT9"|) . 

This completes the third question concerning static 
stability of the bridge. We have found a new catenary 
solution even for bulk charges in the bridge. 



x{t) = t-cf(t), te (0,1 
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with 



d = 2— arcosh- 



and £ to be the solution of the equation 

c = c m (£,6) 
c m (^,b) 
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IV. ANSWER TO QUESTION (IV): 
DYNAMICAL CONSIDERATION 

A. Mass flow of the bridge 

We consider now the actual motion of the liquid in 
the bridge. Here we want to propose the picture that 
possible charges in the water will move according to the 
applied electric field and will drag water particles such 
that a mean mass motion starts. Due to the low Reynolds 
numbers (40-100) for water we can consider the motion 
as laminar and we can neglect the convection term uVu 
in the Navier Stokes equation [25| which reads then for 
the stationary case 



t?V 2 u - Vp + p c E = 0. 



(21) 
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FIG. 3: The upper critical bound for the parameter c accord- 
ing to (I17f> . The inset shows the maximum in dependence on 
the creeping parameter b. 




FIG. 4: The mean mass current through the bridge in depen- 
dence on the electric field and for two different bulk charge 
densities. The thick lines are for a bridge length of 1cm and 
the thin lines for the corresponding length of 2cm. The min- 
imal field strength for stability (|19|l are indicated by corre- 
sponding vertical lines. 



The gradient of the electric pressure © can be given in 
the direction of the bridge by 



e ( £ -l)£ 2 b 
VP = 2L = 2L P9 - 



(22) 



Here we can adopt the stationary pressure since the vis- 
cous pressure is accounted for by the Navier-Stoke equa- 
tion. Assuming that the flow in the bridge has only a 
transverse component which is radial dependent, u(r), 
we can write the Navier Stokes equation (|2"Tj) as 

^(40 +r (^ +c )=° (23) 

with the resulting velocity profile in the direction of the 
bridge 



u(r)-u(R) = 2u \— + < 



r 



(24) 



where R is the radius of the bridge and we have intro- 
duced the characteristic velocity (QJ. Please note that 
we keep the undetermined velocity at the surface of the 
bridge u(R). We will assume in the follwoing that it is 
negligible. The resulting profile (|24|) has the form of a 
Poiseullie flow but with an interplay between forces due 
to bulk charges and dielectric pressure in relation to grav- 
ity. 

The mean current relative to the surface motion is eas- 
ily calculated 



/ = 2irp / drr[u(r) — u(R)] = pvirR 



(25) 



providing the mean velocity of the bridge from (l24l) as 



(26) 



v = u [-+c 



One sees that the ratio of the field-dependent creeping 
height ([2J to the bridge length determines the mean ve- 
locity together with possible dynamical bulk charges de- 
scribed by ([3]). Since we have presently no good control 
over the surface velocity u(R) we approximate it in the 
following as zero. 

Please note that the bulk charge transport described 
by ([3]) leads to Ohmic behavior and the neutral particle 
transport due to dielectric pressure leads to a quadratic 
field dependence condensed in ([2]). The formula (|26|) 
combines the effect of charge transport and neutral par- 
ticle mass transport. It answers the problem raised in 
the introduction how the two pictures can be brought 
together, the one of an incompressible fluids where the 
velocity is dependent of the area and the one of Ohmic 
transport where the velocity is only dependent on the 
electric field. 

The resulting total mass current is given in figure [4] 
The current increases basically with the square of the 
applied field scaled by the bridge length. For additional 
bulk densities the mass flow is higher. 



B. Comparison with the experiment 

To convince the reader about the validity of the veloc- 
ity formula (|26j) we compare now with the mass flow and 
the charge flow measurements. The experimental values 
of Figure 4 in Q are reported to be 40mg/s for a bridge 
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FIG. 5: The necessary applied voltage versus bulk charge 
densities in order to maintain a mass current of 40ml/s. Fol- 
lowing 01 the length of the bridge was L — 1cm and the 
diameter 2.5mm. The result using the flow expression l)26[) of 
the present paper (solid line) is compared to an Ohmic trans- 
port (dashed line). For the latter one the bulk charge has 
been multiplied with 13 orders of magnitude. 



FIG. 6: The necessary applied voltage versus bridge length 
in order to maintain a charge current of 0.5mA. The data are 
from figure 6 of 2]. The result using the flow expression (|26(l 
and a bulk charge of 2.3ng/l (solid line) is compared to an 
Ohmic transport (dashed line). For the Ohmic transport the 
bulk charge has been multiplied with a factor of 3 x 10 3 . The 
same offset of Uq — 8kV is used as in the experiments. 



of lcm length, a diameter of 2.5mm for the stationary 
regime. For this situation we compare in figure [5] the 
results obtained from (|26p with a pure Ohmic transport 
using the lowest-order conductivity expression 



heating easily as 



eN A 



+ cr 



(27) 



where for clean water the conductivity is ctq, A is the 
molecular conductivity of the solved charge (electrolyte), 
and Na the Avogadro constant, see table HI We see that 
our formula (|2^|) leads to a realistic necessary voltage - 
which was 12.5kV in the experiment - even if no bulk 
charge is presented. In contrast, for the Ohmic transport 
one has to assume 13 orders of magnitude higher bulk 
charges to come into the same range. This illustrates the 
advantage of the here presented model. 

Considering the charge transport we do not expect 
such big differences of our model to the pure Ohmic pic- 
ture since the charged particles matters. To this end we 
compare the applied voltage versus bridge length with 
a constant charge current as it was given in figure 6 of 
0. In figure [5] we compare the result from (|26p with the 
pure Ohmic transport. We use a bulk charge of 2.3ng/l. 
In order to obtain a comparable Ohmic result we had to 
multiply the bulk charge with a factor of 3 x 10 3 which il- 
lustrates the difference between our model and the Ohmic 
transport. 

While the difference in charge transport is not very 
significant provided the fact that the conductivity of wa- 
ter varies in the order of 3 magnitudes, the mass flow of 
figure [5] has shown that our result here with (1261) is su- 
perior since it considers the drag of neutral particles due 
to dielectric pressure together with the charge transport. 

Having the current at hand one estimates the Joule 
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(28) 



From figure 5 of Q one sees that the reported increase 
of 10K in 30min would translate into field strengths of 
0.7kV/cm in our calculation. This is much lower than 
our result. We would obtain here 2-3 orders of magni- 
tude higher heating rates. Please note that the cooling 
mechanisms like evaporating and cooling due to water 
flow is beyond the present consideration. Since these are 
probably the major cooling effects in the experiments [26[ 
we cannot compare seriously the theoretical heating rate 
with the experimentally observed ones. 



C. Profile of bridge 

Let us now calculate the profile of the bridge along 
the length. We consider to this end the total mass flow 
of the bridge and neglect the viscous term compared to 
the kinetic energy (which includes part of the convection 
term), uVu = iVu 2 + curlu x u « iVu 2 . Then one 
arrives at the Bernoulli equation 



p^-+pgf{x)+°s{ 



i 



R(x)- 



-p c Ex = p— +a s 



(29) 



Here we have neglected the curvature of the bridge com- 
pared to the curvature due to the radius and have com- 
pared the position-dependent radius R(x) and veloc- 
ity v(x) in the bridge with the situation at the beaker 
(x = 0). The Bernoulli equation (|29p can be rewritten in 
terms of the capillary height (JTJ) and the velocity (J2TJ)) as 



/(*) 



ex 



v 2 (x) 



2.9 



2R(x) 



(30) 



which determines the radius R(x) from the profile of the 
bridge (If 51) and the velocity v(x) if we observe the current 
conservation through an area 



R{x) 2 v{x) = R 2 v. 



(31) 



The results are presented figures [7] and [5] We plot the 
shape of the bridge, the radius and the velocity together 
with a 3D plot. The case of no bulk charges which leads 
to the standard catenary can be found in figure [7] and 
figure [7J shows the situation for extreme bulk charges al- 
most at the stability edge ([20]) . We see a deformation of 
the catenary due to the applied field. This deformation 
is observed, e.g. if an additional field is brought near 
the bridge @, [22|. One sees that the radius is becom- 
ing smaller at one end of the bridge accompanied with 
higher velocities as it is known from falling water pipes 
|27| . The bulk charge leads to deformations of this pro- 
file which are exaggerated in the plot due to the choice 
of unequal scales. 

Interestingly such asymmetry is experimentally ob- 
served [2:] , where after 3 min of operation the asymmetry 
for the bridge of 0.9cm length ranges from a diameter of 
2.1mm to 2.6mm. This is in agreement with the profile 
calculated in figure [8j Also the measured asymmetry in 
the left and right catenary angle Q in glycerine can be 
explained with the present model. 
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D. Dynamical stability 

We turn now to the question of dynamical stability 
of the flow and consider the motion of water together 
with the motion of charged particles characterized by the 
mass mi and charge e^. This charge current is given by 
Ohm's law oE and the corresponding mass current can 
be written 



J, 



-J 



Xi — aE 

Pc 



(32) 



where we introduced the mass ratio of the number of 
charged particles (e.g. NaCl) to the water particles 



#zTO N aCl PclTlt 



(33) 



#u,wih 2 o P^i 
The mass current of the neutral (water) particles are then 

jn = PnVn = (P -Pc)v n = (1 - Xi)pV n (34) 

such that the total mass current reads 

pv = ji + jn = Xi—oE + (1 - Xi)pv n . (35) 

Pc 

The total current (left side) should be larger than the 
current only from the charged particles (last term on the 
right side). However the velocity of charged particles, 
aEj p c should be larger than the velocity of the dragged 




FIG. 7: The center of mass coordinate (above), the radius 
(middle) and the velocity (bottom) together with the 3D plot 
of water bridge (in cm) for no bulk charges c = 0. The pa- 
rameter are b = 1.5cm and according to table [I] Please note 
the different length scales in x and y, z direction. 



water molecules v n and therefore larger than the mean 
velocity v of the mass motion. Together with (|2"oT) this is 
expressed by the inequalities 



oE b \ aE 

>Uq[—+C) > Xi 

p c \L J p c 



(36) 
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FIG. 8: The center of mass coordinate (above), the radius 
(middle) and the velocity (bottom) together with the 3D plot 
of water bridge (in cm) with bulk charges c = 1. The param- 
eter are b = 1cm and according to table [I] 



which gives an upper and lower bound on the possible 
mass motion created by the drag of particles due to the 
force on charged particles. 

If we now take into account the dependence of the con- 
ductivity on the density of the solved ions in water we 
can find a condition on possible bulk charges in water to 
maintain a stable bridge. To this aim we consider very 
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FIG. 9: The range of possible water bridges for an electric 
field of E = 0.64kV/cm. The upper limit is due to the static 
stability condition ()20p and the lower cut is due to the dy- 
namical condition (|37[l . The bulk-charge-free condition is the 
upper straight line. 



small charge densities solved in water which allows to 
consider the lowest order dependence of the conductivity 
on the bulk charge concentration (j2"T)l . 

Noting the charge-density dependencies of Xi , b and c 
via ([231) , © an d © one obtains from ([55)1 the dynamical 
restriction on possible bulk charges 



p c e pi- p 2 ±J{pi- P2) 2 +pj 



p c (l- 2p 2 /pi) > pl/pi - 2pi 



(37) 



with the auxiliary densities 
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(38) 



The results for NaCl in water (table [TJ are plotted in 
figures CTTTJl The static stability condition (|19p gives the 
upper and charge-density- independent limit in figure [TO] 
The static condition (|2TJ)) with bulk charges leads to the 
border of maximal densities on the right side which agrees 
with (|19[) at zero densities, of course. The lower minimal 
length of the bridge at a given field strength and bulk 
charge is provided by the dynamical condition (|37p . For 
no bulk charge the possible range of lengths of the bridge 
starts at zero and is limited by the upper length (|19p . 
If there are charges present, there is a minimal length 
required to have a stable bridge. 

From the 3D plot in figure [TU] one can see that for 
finite charges and for fixed bridge lengths there is a lower 
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FIG. 10: The range of possible water bridges in dependence 
on the bridge length, the electric field and the electrolyte bulk 
charges. 



and an upper critical field where bridges can only be 
stable. From the experiments Q it is seen that the bridge 
forms jets for fields higher than 15kV/cm and therefore 
becomes unstable. With a bridge length of 0.5cm this 
translates into a bulk charge of 4ng/l according to our 
found boundary conditions. This is in agreement with 
the value needed to reproduce the flow measurements 
described in section IV. B. 



V. SUMMARY 

The formation of water bridges between two vessels 
when an electric field is applied has been investigated 
macroscopically. Electrohydrodynamics is sufficient to 
describe the phenomena in agreement with the experi- 
mental data. The four necessary parameters which are 
build up from microscopic properties of the charged liq- 
uid are the capillary height ([T|), the creeping height ((2J, 
the dimensionless ratio between field and gravitational 
force density ([3]), and the characteristic velocity (U). 

As new contribution to the discussion, an exact solu- 
tion has been found of a charged catenary. This leads 
to a static stability criterion for possible charges in the 
liquid dependent on the applied field strengths and on 
the length of the bridge. With no bulk charges present 
the maximal bridge length is determined and no minimal 
length occurs. This changes if bulk charges are present. 
Then also a minimal length is required. However, only 
very small concentrations of bulk charges are possible and 
the bridge is easily destroyed when bulk charges exceed 



50 ng/1. As a further an asymmetric profile in the diam- 
eter along the bridge is obtained which was observed by 
asymmetric heating. 

For the dynamical consideration a picture is proposed 
of dragged liquid particles due to the motion of the 
charged ones besides the ponderomotoric forces due to 
the dielectric character of the liquid. The resulting con- 
sideration of dynamical stability restricts the possible pa- 
rameter range of bridge formation. The resulting mass 
flow combines the charge transport and the neutral mass 
flow dragged by dielectric pressure and is in agreement 
with the experimental data. 

The presented simple classical theory applies for 
charged liquids as long as the Reynolds number is such 
low that laminar flow can be assumed. 
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Appendix A: Solution of charged catenary 

Here the drivation of the charged catenary [24| is 
shortly sketched. We solve the variation problem (IT41 



1^ 



(Al) 



with the functional J-{x) — pg [f(x) + b — cx\ \J\ + f'(x) 2 
and the boundary conditions /(0) = f(L) = 0. 
It is useful to introduce 



t(x) = f(x) +b — ex 



such that 



F{x) = pgt{x)jl + [t'(x) + > 
The corresponding Lagrange equation 
d ST ST 



dx St'(x) St{x) 
possesses a first integral 
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T = const 
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(A2) 
(A3) 

(A4) 
(A5) 



where we introduced the first integration constant £ in a 
convenient way. 
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The resulting differential equation 



t{x)[ct'(x) + 1] = ^t'{x) 2 + (ct'(x) + l) 2 (A6) 
with x = x(l + c 2 ) is solved in an implicit way 



t(x) = £cosbJ - 



x + ct(x) — cb + —d 



(A7) 



with a second integration constants d. The profile is 
therefore given by the implicit equation 



f(x) = ex — b + £ cosh < - 



x + cf(x) + —d 



(A8) 



The boundary condition /(0) = leads to the deter- 
mination of the integration constant 



£ fb 

d = 2— arcosh I - 



(A9) 



in terms of the yet unknown £ constant. The solution 
can be written with the help of (|A9|) as 



f(x) — cx + t; <cosh 



x + cf(x) Ld 

I 2e 



-cosh I — - 
\2€ 



(A10) 



The boundary condition f(L) = lead to the determi- 
nation of the remaining constant £ to be the solution of 
the equation 



■ — sinh — 
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b ■ u L 
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2£ U 2 £ 
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Finally we can rewrite the implicit solution (|A10j) in 
parametric form. Therefore we choose as parameter 
t = x + cf(x) which runs obviously through the inter- 
val t £ (0, L) and we obtain the solution (fT5|) 



/(*) 



1 



1 + c 2 
x(t) = t-cf(t), 



ct+C 



cosh ( ) — cosh I — ) 



te (0,L). 
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